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1.  Introduction. 

Physical  effort  at  high  altitude  can  result  in  serious  complications  in  the  human  respiratory  system  even  for  healthy 
and  well-trained  persons.  For  people  with  a pulmonary  disorder,  already  a stay  at  moderate  altitude,  or  transportation 
by  air  (e.g.  in  case  of  medical  evacuation)  can  lead  to  significant  problems  caused  by  hypoxemia. 

In  both  circumstances  the  oxygen  tension  of  the  inspired  air  (PI,02)  and  therefore  of  the  arterial  blood  (Pa,02)  drops 
substantially.  For  healthy  people  - if  at  rest  - this  causes  no  considerable  problems  because  their  arterial  oxygen 
saturation  at  such  altitudes  still  will  be  at  the  horizontal  part  of  the  oxyhaemoglobin  dissociation  curve. 

Problem  area 

Patients  with  COPD  and  interstitial  lung  diseases  have  mild  to  severe  obstructive  and/or  restrictive  impairment 
of  the  lung  function.  So  they  may  deal  with  unequal  ventilation,  ventilation-perfusion  mismatch,  diffusion  disorders, 
ventilatory  restriction,  and  changes  in  arterial  blood  gas  tensions.  In  this  kind  of  lung  diseases  there  is  a diminished 
reserve  of  gas  exchanging  surface  because  of  degeneration  of  alveoli  and/or  fibrosis  of  the  lungs. 

For  the  individual  patient  with  COPD  or  interstitial  lung  diseases  the  Pa,02  at  altitude  can  not  be  predicted 
from  the  sea  level  value  (Ref.  1).  Even  with  normal  blood  gas  tensions  at  sea  level  severe  drop  in  oxygen  saturation 
may  occur  at  altitude.  Predicting  factors  seem  to  be  a FEV1  (forced  expiratory  volume  1st  second)  < 30  % VC  (vital 
capacity)  and/or  a significant  decrease  of  the  Pa,02  with  comparatively  low  exercise.  Other  empirical  formulas  to 
predict  Pa,02  at  altitude  from  sea  level  conditions  have  been  derived  by  Gong  (Ref.  2)  and  Dillard  (Ref.  3).  These 
statistical  formulas  indicate  more  or  less  average  responses  but  are  not  suitable  for  the  predictions  for  the  individual 
patient. 

HAST  (high  altitude  simulation  test)  may  provide  more  clearness  in  the  individual  for  the  extent  of  the  gas 
exchanging  surface.  An  important  difference  of  this  test  with  regard  to  reality  is  the  constant  barometric  pressure  and 
the  comparatively  limited  duration  (30  min).  Besides,  this  test  is  not  without  any  risk  for  the  patient  (hypoxemia 
induced  arrhythmias). 

Approach 

In  this  paper  we  explore  the  modeling  of  the  human  lung  gas-exchange  under  hypobaric  conditions.  At  first 
the  concept  of  the  ventilation-perfusion  mismatch  is  applied  to  model  the  pulmonary  performance  of  a healthy  person  at 
rest  at  increasing  altitudes.  Model  results  are  compared  with  experimental  data  from  literature  at  varying  altitude.  The 
model  is  also  used  for  a patient  with  COPD  at  moderate  altitude.  The  results  illustrate  the  limited  applicability  of  the 
mismatch  model.  Inclusion  of  diffusion  effects  improves  the  results.  Only  a preliminary  validation  of  the  model  has 
been  attempted  yet.  The  focus  is  predominantly  on  understanding  and  modeling  of  possible  underlying  physics. 

2.  The  ventilation-perfusion  mismatch  model. 

Compartments 

The  algorithm  for  the  exchange  of  02  and  C02  is  based  on  the  subdivision  of  the  human  lung  into  a number  of 
compartments.  Several  units,  consisting  of  alveoli  and  the  network  of  capillaries,  are  lumped  together  in  each 
compartment.  One  assumes  that  the  conditions  for  all  units  in  one  compartment  are  the  same.  The  ventilation-perfusion 
ratio  vVq’  for  a compartment  dictates  the  transfer  of  02  and  C02  as  a function  of  their  mixed  venous  values  PV,02, 
PV,C02  and  the  oxygen  inspired  presure  PI,02.  Two  examples  for  a range  of  values  of  vVq*  are  shown  in  figure  1. 
Curve  A shows  the  variation  of  end-capillary  values  of  partial  oxygen  pressure  Pec,02  versus  Pec,C02  for  sea  level 
atmospheric  conditions,  curve  B for  hypobaric  conditions.  Each  point  on  a curve  represents  a compartment  with  a 
certain  ventilation-perfusion  ratio. 

The  basis  for  the  construction  of  these  curves  is  that  complete  equilibration  of  02  and  C02  has  taken  place 
between  the  capillary  blood  flow  and  the  alveoli.  For  C02  this  is  not  an  issue  because  its  transfer  through  the  alveolar 
membrane  takes  place  very  rapidly.  However  for  02  this  is  only  true  if:  a)  the  blood  barrier  resistance  of  the  membrane 
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is  not  too  high,  b)  the  contact  time  of  the  capillary  flow  with  the  alveolar  tissue,  is  not  too  small,  and  c)  the  initial  PO2 
difference  over  the  membrane  is  not  small  as  might  occur  in  hypobaric  conditions.  We  will  discuss  items  b)  and  c)  later 


The  method  of  Rahn  and  Riley  & Cournand  described  by  Farhi  (Ref.  4)  is  used  to  find  the  points  of  figure  1. 
The  equations  are  given  for  completeness: 


k 


02  = v'n *(P/, 02  - PA, 02)*2-  = q\  02  = q\ *{cec, 02  - cV, 02 ) 


n =1  ....N 


(1) 


It  states  that  for  each  compartment  n the  amount  of  loss  of  O2  in  the  gaseous  state,  v’n02,  equals  the  transfer  of 
02  to  the  capillaries,  q’n02.  The  concentration  in  the  blood  at  the  end  of  the  capillary  is  indicated  by  cec,02,  at  the 
venous  point  cV,02 ; k is  a conversion  factor,  v’n  and  q’n  are  the  ventilation  and  perfusion  of  compartment  n. 

The  uptake  of  CO2  into  the  alveolar  gas  from  the  blood  is  expressed  as: 


v\  C02  = v\  *PA,  C02*-  = q\  C02  = q'„  *(cec,  C02  - cV,  C02 ) (2) 

k 

The  equations  for  the  concentrations  c02  and  cC02  are  taken  from  West  and  Wagner  (Ref.  5).  They  are  both  non-linear 
functions  of  Hb,  hematocrit,  P02,  PC02  and  body  temperature.  By  the  assumption  of  fully  equilibration  the  end- 
capillary values  for  P02  and  PC02  are  equal  to  their  corresponding  alveolar  values: 


PA,02  = Pnec,02 
PnA,C02  =P„ec,C02 


Ventilation-perfusion  distributions 

In  analogy  to  the  model  of  West  and  Wagner  we  assume  a discrete  lognormal  distribution  of  ventilation  and 
perfusion  over  the  compartments: 


Q *A|  * 
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(4) 

(5) 


V’A  and  Q’  are  the  total  alveolar  ventilation  and  pulmonary  perfusion.  In  figure  2 an  example  is  shown.  The  v’n  and  the 
q’n  are  plotted  versus  their  logarithmic  ratio  = log  (v’n  / q’n)-  The  discrete  compartments  are  centered  at  %n.  All 
compartments  are  equidistant  in  ^ and  have  equal  widths  A^.  The  parameter  a functions  in  these  formulae  in  two  ways. 
In  both  the  single  distributions  for  v’n  and  q’n  it  is  a measure  of  the  deviation  from  the  means  juv  and  pq.  Its  square  a2 
also  separates  both  peaks.  From  (4)  and  (5)  it  follows: 

V,  =log(9-)+f  (6) 

M,=log(^)-ic7! 

We  will  use  G as  the  parameter  representing  the  ventilation-perfusion  mismatch.  For  a healthy  young  person  the 
mismatch  a is  small  Then  both  peaks  for  ventilation  and  perfusion  lie  close  together  and  the  width  a of  their  bell 
shaped  distributions  is  also  small.  We  assume  in  case  of  COPD  that  the  lung  degradation  can  be  expressed  in  an 
increase  of  g.  An  example  of  a pure  variation  in  mismatch  is  shown  in  figure  3.  The  growth  of  g will  lead  to  a steady 
decrease  of  Pa,02  and,  to  a lesser  degree,  to  a steady  increase  of  Pa,C02.  The  responses  of  the  pulmonary  ventilation  to 
an  increase  of  Pa,C02  or  a decrease  of  Pa,02  are  not  taken  into  account  yet.  The  metabolic  indices  V,02  and  V,C02  are 
kept  constant.  Therefore  the  values  of  alveolar  partial  pressures  of  02  and  C02,  averaged  over  all  compartments,  are 
constant  as  well.  By  the  assumption  of  fully  equilibration  one  can  refrain  from  the  inclusion  of  alveolar  membrane 
properties  and  from  the  effect  of  reaction  speed  of  the  02  molecules  with  Hb. 


Circulation 

We  use  two  extra  equations  due  to  Fick’s  principle.  They  describe  the  consumption  of  02  and  the  production 
of  C02  at  the  systemic  side  of  the  circulation: 


V'02  = CO  ■ (ca,  02-cV,02) 
V'C02  = CO  ■ (cV,  C02  - ca,  C02 ) 


(7) 
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The  first  equation  states  that  the  total  02  consumption  equals  the  arterial  and  mixed  venous  concentration  difference 
times  the  cardiac  output  CO.  The  second  equation  expresses  the  equivalent  for  C02  production.  By  the  summation  of 
equations  (1)  and  (2)  over  all  compartments  including  the  right-to-left  shunt,  the  same  quantities  have  to  found  for  the 
02  and  C02  transfer  in  the  total  lung. 

Baseline  conditions 

To  estimate  the  mismatch  parameter  a for  a person,  our  point  of  departure  is  the  knowledge  of  (measured)  arterial 
partial  gas  pressures  Pa,02  and  Pa,C02.  The  cardiac  output  CO  under  resting  conditions  is  estimated  from  the  body 
surface  area  BSA.  Also  the  pulmonary  shunt  fraction  psh  has  to  be  measured  or  estimated.  The  metabolic  rates  V’02 
and  V’C02  are  usually  measured  as  a part  of  exercise  testing.  For  the  perfusion  distribution  to  be  fully  known  only  o 
has  to  be  determined.  The  ventilation  distribution  requires  also  o but  in  addition  the  total  alveolar  ventilation  V’A.  V’A 
will  be  determined  simultaneously  with  a.  The  non-linear  system  of  equations  (1)  to  (7)  is  solved  by  Broydn’s 
procedure  (Ref.  6)  for  chosen  values  of  the  interval  and  the  total  number  of  compartments  N.  The  Broydn  procedure 

is  a fast  and  robust  method  that  finds  the  roots  in  analogy  with  the  secant  method  working  in  two  dimensions.  Hence  it 
does  not  require  the  elaborate  evaluation  of  a Jacobian  matrix.  The  matching  procedure  delivers  the  four  quantities: 
mismatch  parameter  a,  alveolar  ventilation  V’A  and  the  mixed  venous  partial  pressures  PV,02  and  PV,C02. 

To  display  the  interactions  of  the  solution  process  consider  figure  4.  It  shows  Pa,02  and  Pa,C02  as  functions  of  o 
and  alveolar  ventilation  V’A.  The  mixed  venous  partial  pressures  are  not  visualised.  According  to  physical  feeling 
Pa,C02  is  a decreasing  dominant  function  of  V’A  and  it  is  a weakly  increasing  function  of  a.  Pa,02  depends  equally  on 
both  variables.  In  general  it  increases  with  V’A  and  it  decreases  with  a.  For  a measured  pair  of  arterial  partial  pressures 
of  02  and  C02  the  corresponding  values  of  V’A  and  a can  be  found  by  the  Broydn  procedure.  From  a first  guess  of 
V’A  and  a with  probably  incorrect  corresponding  values  of  Pa,02  and  Pa,C02  a trajectory  in  the  V’A-g  plane  is 
traversed  that  tends  to  approach  the  measured  values  of  Pa,02  and  Pa,C02. 

Hypobaric  conditions 

In  our  model  no  explicit  reference  is  made  to  the  atmospheric  pressure,  hence  we  study  the  consequences  of 
hypobaric  conditions  by  changing  the  02  fraction  Fl,02  of  the  inspiratory  gas.  As  a starting  point  we  assume  that  the 
mismatch  parameter  o and  cardiac  output  and  the  right-left  shunt  are  kept  the  same  and  hence  the  perfusion  distribution 
is  kept  the  same.  However  the  ventilation  is  allowed  to  increase  as  a response  to  hypoxia: 

V'E  = V'b  E - * ( Sba , 02  - Sa,  02  )p  (8) 

Where  V’E  is  the  pulmonary  minute  ventilation,  superscript  b indicates  quantities  measured  or  calculated  under  basic 
(i.e.  sea  level  atmospheric  and  resting)  conditions.  Gsa  is  the  sensitivity  (a  negative  number)  of  the  minute  ventilation  to 
arterial  oxygen  saturation.  For  an  exponent  p = 1 Gsa  is  in  the  range  of  -1.55  ±0.98  (Ref.  7).  A hypocapnic  effect  on  GSa 
is  not  taken  explicitly  into  account.  The  real  response  to  hypoxic  conditions  depends  also  on  the  speed  by  which 
hypoxic  conditions  are  imposed  and  on  the  additional  control  of  alveolar  PA,C02.  (Ref.  8). 

From  the  minute  ventilation  the  alveolar  ventilation  can  be  found  by: 

V'A  =V'E  - f(V'E  )*VD  (9) 

The  breathing  frequency  f is  a weak  function  of  V’E  itself;  VD  is  the  physiological  dead  space.  The  system  of  equations 
(1)  to  (9)  can  be  solved  for  the  venous  and  arterial  blood  gasses  using  Broydn ’s  procedure.  The  mismatch  a is  kept 
fixed  now,  but  the  ventilation  V’A  is  allowed  to  adapt  itself  to  hypoxic  conditions  according  to  (8)  and  (9). 

3.  Application  of  the  model. 

Healthy  person 

At  first  the  mismatch  model  is  applied  to  a healthy  subject  starting  at  sea  level  and  decreasing  its  FI,02  to  an 
equivalent  altitude  of  9 km  (Figure  5).  From  the  barometric  pressure  for  standard  atmosphere  the  P1,02  is  found  and 
used  as  an  input  of  the  model.  A low  value  of  mismatch  a is  assumed  that  is  kept  fixed.  The  pulmonary  ventilation  was 
allowed  to  respond  to  hypoxia  according  to  (8)  with  an  exponent  Vi.  This  weak  response  is  a result  of  hypoxia,  damped 
by  the  decreased  PA,C02.  The  results  for  arterial  pressures  are  quite  comparable  with  actual  measurements  of  a 
simulated  ascent  of  the  Everest  (Ref.  9).  In  the  measured  realistic  conditions  the  alveolar-arterial  P02  difference 
increases  with  altitude.  At  9 km  this  amounts  to  a P(A-a),02  = 6 mmHg.  This  might  be  attributed  to  the  lack  of  complete 
equilibration.  The  results  of  the  present  model  show  a steady  decrease  of  the  alveolar-arterial  difference,  practically 
approaching  zero  for  the  higher  altitudes.  This  is  due  to  the  pattern  of  the  02  dissociation  curve.  At  low  02  saturation 
the  same  shunt  gives  a smaller  effect  on  Pa,02  versus  PA,02. 
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Patient  with  COPD 

So  far  things  seem  to  compare  reasonably  well.  However  application  of  the  present  method  to  the  data  of  a 
COPD  patient,  put  forward  by  the  St.  Antonins  Hospital,  showed  results  deviating  from  the  measured  data  (Fig.  6).  It 
seems  that  the  ventilation-perfusion  mismatch  is  not  the  only  source  attributing  to  the  failure  of  this  subject  to  cope 
with  the  lowered  PI,02.  One  of  the  possibilities  is  discussed  below. 


Diffusion  limitations 

Many  patients  suffering  severe  COPD  have  lost  a great  deal  of  their  oxygen  transfer  capabilities.  Under 
normal  resting  conditions  this  does  not  necessarily  lead  to  hypoxemia.  But  in  exercise  and  under  hypobaric  conditions 
this  can  be  a cause  of  hypoxemia.  The  transfer  properties  of  the  lung  are  evaluated  clinically  by  measuring  a.o.  TL,CO 
that  is  defined  as: 


TL,CO 


V'CO 
PA, CO 


(10) 


Which  measures  the  uptake  of  carbon  monoxide,  CO.  See  e.g  reference  10  for  a careful  analysis  of  TL,CO 
measurements.  CO  has  such  a large  affinity  for  Hb  that  its  capillary  pressure  can  be  assumed  zero  and  PA, CO  can  be 
considered  as  the  driving  pressure.  Therefore  TL,CO  characterizes  the  properties  of  the  blood  gas  barrier.  It  is 
proportional  to  the  membrane  surface  area  and  the  solubility  of  CO  and  is  inversely  proportional  to  membrane 
thickness  and  the  square  root  of  the  molecular  weight  of  CO.  A reduction  of  membrane  area  by  COPD  will  lead  to  a 
proportional  reduction  in  transit  time  of  the  blood  in  the  capillary  bed  if  the  same  cardiac  output  has  to  be  maintained. 
All  other  membrane  properties  are  assumed  the  same.  A decrease  of  transit  time  tc  is  no  problem  if  the  oxygen  equil- 
ibration time  te  remains  sufficient  small.  However  for  lower  values  of  PA, 02  this  is  not  possible.  In  figure  7 three 
examples  are  shown  of  alveolar  capillary  P02  differences  during  the  passage  of  blood  along  alveolar  tissue.  The  curves 
are  obtained  from  approximations  of  the  diffusing  equation  for  oxygen: 


—7 ± = g*D02*(PA,02-Pc,02 ) <n> 

at 

In  which  t is  the  time  coordinate  for  a blood  parcel  moving  along  a capillary,  g is  a geometrical  shape  parameter  and 
D02  is  the  02  membrane  permeability  factor.  The  02  concentration  c02  in  the  capillary  vascular  system  depends  on  the 
local  gas  pressures  Pc,C02  and  Pc,02.  All  parameters  are  taken  for  healthy  non-degraded  tissue.  Normal  value  for  the 
contact  time  is  tc  = 0.76  sec  (Ref.  11).  Equilibration  time  te  is  small  for  a large  alveolar- venous  oxygen  pressure 
difference  and  a PA,02  on  the  horizontal  part  of  the  02  dissociation  function  (upper  curve),  te  increases  strongly  with 
decreasing  alveolar  oxygen  pressures  (curves  PA,O2=90,  60).  For  a fixed  tc  < te , the  alveolar-end-capillary  difference 
P(A-ec),02  is  small  for  high  PA,02  because  te  is  small.  P(A-ec),02  increases  for  intermediate  values  (e.g.  PA,02  = 90) 
because  te  is  increasing.  However  P(A-ec),02  decreases  again  for  smaller  PA,02,  because  the  alveolar-venous  difference 
is  decreasing.  But  it  will  not  vanish.  For  lower  PA,02  there  is  always  a condition  where  no  equilibration  can  occur 
because  te  will  be  much  larger  than  any  tc. 

In  our  system  of  equations  the  first  of  equation  (3)  is  replaced  by: 


Pn  ec,  02  (tc)  = PnA,02-Pn(A-  ec),  02  [g,  DQ2 ; tc]  n=1  N (12) 

Where  the  last  term  follows  from  an  approximate  solution  of  (1 1)  for  each  compartment  and  for  the  yet  unknown  initial 
condition  PV,02  and  the  compartment  dependent  PA,02.  Broydn’s  procedure  again  is  applied  but  now  with  N additional 
non-linear  equations  (12). 


Patient  with  COPD  and  diffusion  limitations 

Baseline  results  for  FI,02  =0.21  of  our  COPD  patient  are  shown  in  figure  8 A.  It  displays  alveolar  and  end-capillary 
values  of  P02  versus  the  ventilation-perfusion  ratio.  From  Tl,C0  data  it  is  inferred  that  the  average  contact  time  was 
reduced  to  about  tc=03&  sec.  For  alveolar  PA,02  = 120  mmHg  and  higher  the  equilibration  takes  place  before  the  end 
of  the  capillary  is  reached.  For  lower  alveolar  values  this  is  not  the  case  and  alveolar-end-capillary  differences  are 
shown.  This  has  however  not  a strong  detrimental  influence  because  the  ventilation  as  well  as  the  perfusion  are  smaller 
in  that  region.  Figure  8B  shows  the  distributions  for  FI,02  = 0.15.  Over  the  whole  range  of  ventilation-perfusion  ratios  a 
substantial  alveolar-end-capillary  difference  occurs  in  P02.  The  smaller  differences  now  appear  at  the  lower  side  of  the 
PA,02  because  here  the  smaller  alveolar-venous  difference  starts  to  reduce  the  alveolar-end-capillary  differences. 
Figure  9 depicts  the  course  of  the  blood  gas  values  starting  from  normal  FI,02  = 0.21  to  gradual  lower  FI,02.  A 
reasonable  match  with  HAST  measurements  is  obtained. 


4.  Concluding  remarks. 

At  hand  of  a healthy  person  and  one  with  COPD  the  application  of  the  ventilation-perfusion  mismatch  model 
of  the  human  lung  has  been  investigated.  The  model  can  not  predict  the  pulmonary  performance  at  moderate  altitude 
for  a COPD  patient.  It  is  likely  from  diffusion  data  that  the  reduced  contact  time  of  the  blood  flow  in  the  alveolar  bed 
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may  cause  the  remaining  discrepancies.  To  include  such  possible  effects,  the  model  has  been  extended  with  a capillary 
O2  equation  that  accounts  for  alveolar-end-capillary  P02  differences.  A better  match  with  the  Pa,02  drop  observed  in 
HAST  data  appears  to  justify  these  model  extensions.  The  validation  of  the  improved  model  requires  however  the 
application  to  a greater  number  of  subjects.  This  will  be  pursued  in  future  work.  A satisfactory  prediction  of  oxygen 
saturation  at  altitude  has  to  take  into  account  more  person  specific  characteristics.  Hence  the  extension  of  the  model 
might  not  be  limited  to  ventilation-perfusion  mismatch  and  diffusion  effects  alone. 
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Figure  1.  Alveolar  PC02  versus  P02  for  a subject  at  sea  level  (curve  A)  and  at  h = 2500  m (curve  B).  Each 
point  on  a curve  represents  a compartment  with  a certain  v’/q’  combination.  At  the  venous  points  (V) 
v7q’=0  and  at  the  inspired  points  (I)  v7q’  approaches  infinite. 


Figure  2.  Ventilation  v (L/min)  and  perfusion  q’  (L/rmn)  lognormal  distributions  as  function  of  their  ratio 
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Figure  3.  Effect  of  mismatch  parameter  a oil  alveolar  (PA,02,  PA,C02)  and  arterial  partial  pressures  (Pa,02, 
Pa,C02). 


Figure  4.  Lines  of  constant  Pa,02  (thin  lines  with  grey  fill)  and  of  constant  Pa,C02  (dotted  heavy  lines)  in 
mmHg  as  function  of  mismatch  parameter  a and  alveolar  ventilation  V’ A.  For  a pair  of  measured  values  of 
Pa,02  and  Pa,C02  the  characteristics  V’A  and  a can  be  found. 
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Figure  9.  Effect  of  lowering  Fl,02  on  Sa,02,  Pa,02,  Pa,C02  and  V’E.  At  the  bottom  altitude  is  given  in  km. 
Open  symbols  represent  HAST  measurements  of  Sa,02?  Pa,02,  Pa,C02  of  COPD  patient. 

Crosses  are  Pa,02  from  statistical  formulae  of  Gong  and  Dillard  (Refs.  2,3). 

Lines  are  computed  by  ventilation-perfusion  mismatch  method  including  non-equilibration  of  capillary  P02. 
Initial  conditions  at  Fl,O2=0.21  are  derived  from  measurements. 
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